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FOREWORD 


This annual technical report documents research performed from 1 February 1973 
tc 31 January 197^ under Office of Naval Research Contract N00014-73-C-0281. The 
research program is beirg conducted at United Aircraft Research laboratories (UARL), 
East Hartford, Connecticut O 0 IO 8 . Mr. C-overt Flohil is serving as the ONR Scientific 
Officer. 

This report is issued as UARL Report N9U620-2. 
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Development of Optimal Control Modes 
for Advanced Technology Propulsion Systems 

SUM-iARY 


A nonlinear multivariable feedback controller was defined for the idle to 
military operating regime (9 to 100 percent thrust) of the Pratt & Whitney Aircraft 
F401 variable cycle turbofan engine. The analytical design involved (1) linearizing 
the F401 engine dynamics about five steady-state operating points between idle and 
military thrust, (2) applying linear optimal control synthesis methods at each point 
and (3) combining the five optimal linear controllers into a single nonlinear con¬ 
troller which has feedback gains that are scheduled with high compressor speed. 
Variable fan, compressor and exhaust geometries as well as main burner fuel flow are 
coordinated by the controller to achieve rapid engine dynamic response. 

Engine accelerations from idle to military thrust levels with the defined non¬ 
linear (optimal) controller were computed using a detailed digital nonlinear simula¬ 
tion of the engine. These accelerations were compared with those obtained using a 
conventional controller designed to provide rapid thrust response. The optimal con¬ 
troller provided significantly faster F401 engine model acceleration from idle to 
military thrust, without exceeding temperature or stability margin constraints. 
Moreover, the optimal controller moved exhaust and fan geometries in a significantly 
different manner than t-ne conventional controller. 

An integrated inlet-engine controller was also defined for the F401 engine and 
a mathematical model of an internal compression supersonic inlet. The inlet had 
variable throat and bypass geometries and wpj interface., with the F401 engine simu¬ 
lation for a flight condition of 40,000 ft and Macb ?,2. Closed-loop propulsion 
system respc.xoe to simulated afterburner ignitior as -'■vaiuated for separate inlet- 
engine controls and for the integrated controller. 

The integrated controller provided closed-loop regula'. on xhich w is as good 
as or better than the separate controls for all critical inlet and engine variables. 
Also, the improved regulation was accomplished using consi ^erably sr uier cos ;rol 
positions and rates. The improvement in inlet-engine dynamic performance with the 
integrated controller was due primarily to cross-feedback from the engine state 
variables to the inlet control variables. 
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RESULTS AND CONCLUSIONS 


1. A nnlinear multivariable feedback controller was developed for operation of 
the Fuel variable geometry turbofan engine between idle and military thrust condi- 
tiens '9 percent and IX percent thrust, respectively). The piecewise-linear/:iece- 
wise-::timal approach to nonlinear control synthesis, developed and reported under 
a ; revi us CNR contract (Ref. 1’, vac used to define the controller. Thrust restonse 
for nonlinear multivariable control was approximately 30 percent faster than that 
for tr.e conventional ccntrcller, i.e., the 98-percent thrust point was reached in 
2.32 sec versus 3*^0 sec for conventional control. Beth controllers avoided engine 
overtemperatures and provided similar fan and compressor stability margins. 


2. The nonlinear multivariable controller changed fan inlet guide vane positions 
and Jet exhaust area in a manner significantly different from the conventional con¬ 
troller. For the optimal control mode, fan inlet guide vane changes led the steady- 
cta‘r .r ; r 'i.*,.. fc \e rule vane versus fan r ec; - she bile programmed for rt^n-lv- 
state F^Ol operation), and for conventional control they lagged the steady-state 
schedule. Jet exhaust area, with the optimal controller, did not move far from its 
initial position until the fan stability margin had begun to increase. At that time, 
exhaust area decreased rapidly tc its bottom position limit and then increased rapid- 






its final steady-state value. In the conventional control mode, exhaust area 


first increased rapiily, remained constant for a time, and then decreased to its 
steady-state value. 


3. A mathematical model of an internal compression supersonic inlet was interfaced 
with the FUoi engine simulation at a flight condition of U0,000 ft and Mach 2.2, 
and an integrated inlet-engine controller was designed for that operating point. 
Integrated control also was synthesized using the piecewise-linear/piecewise-optimal 
techniques. The integrated controller regulated all critical inlet and engine 
variables as well as, or better than, the separately designed system. The improved 
performance was accomplished using considerably smaller control positions and rates, 
and was due to inlet-engine cross-feedback in the controller. 
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HfEROKJCTION 


Advanced technology propulsion systems consist of inlet-engine combinations 
which incorporate a variety of variable geometry features. Traditional control 
synthesis procedures appear to be inadequate to exploit the full performance capa^- 
bilities cf advanced technology propulsion systems because the traditional tech¬ 
niques do not properly treat the parameter interaction which exists in these non¬ 
linear multivariable systems. However, modern methods based on a state variable 
description of system dynamics and the use of optimization procedures to generate 
feedback control provide the necessary analytical techniques. 

The piecewise-linear/piecewise-optimal control synthesis approach, in partic¬ 
ular, is attractive because it provides a means to analyze the nonlinear aspects of 
the propulsion system controller design problem. In addition, there are no compro¬ 
mises in its ability to handle the multivariable properties of the system. This 
control design concept was developed by UARL in a previous program with ONR (Ref. 1). 
Optimal multivariable feedback controllers are developed for each point about which 
the dynamic system to be controlled is linearized. N nlincar feedback • \tr 1 is 
then constructed by scheduling the optimal linear controllers with system state. 

The nonlinear characteristics of the dynamics influence the analysis only in the 
initial phase during piecewise linearization and in the final phase when the non¬ 
linear controller is implemented and evaluated. In a previous program (Ref. 2), 
this technique was applied to a detailed nonlinear model of the Pratt & Whitney 
Aircraft f 401 variable cycle turbofan engine and an optimal multivariable feedback 
controller was defined for high power engine operation (80 to 100 percent thrust). 
Calculated F401 engine response using this controller was significantly faster than 
the response for a state-of-the-art conventional control*. 

The principal objective of the studies described in this report was to define 
and evaluate an optimal multivariable controller capable of accelerating the F401 
engine model from idle (9 percent thrust) to military (100 percent thrust) condi¬ 
tions. Another maj^r objective was to determine performance advantages and char- 
acterist_cs of an inlet-engine feedback controller designed by treating the inlet 
and engine as a single integrated dynamic system. 


*The conventional control used for comparison in Ref. 2, as well as in this report, 
is the P&WA fast accel controller. 
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SYNTHESIS OF AN OPTIMAL FEEDBACK COI7TROL SYSTEM 


The application to the F401 engine of the piecewise-linear/piecewise-optimal 
method for synthesizing multivariable feedback control is described in this section. 
The objective was to synthesize a multivariable controller which would provide 
significantly faster accelerations than those resulting for a conventional control 
mode, while satisfying temperature, stability, control position aid control actua¬ 
tion rate constraints. Engine input, state and output variables which are appro¬ 
priate for the synthesis procedure are introduced initially, and the linearization 
of engine dynamics is described. Linear optimal control methods used to compute 
feedback gains are outlined next, and the development of the nonlinear feedback 
controller from the sets of linear gains is discussed. Finally, the particular 
optimal multivariable controller developed for the F401 engine is detailed, and 
calculated engine response for the optimal multivariable controller is compared 
with that which can be obtained for a conventional controller. 


Multivariable Engine Dynamics 

The first step in applying the piecewise-lineer/piecewise-optimal control 
technique is to define engine control variables, state variables, and output 
variables. The FhOl engine has variable exhaust, fan, and compressor geometries. 
Accordingly, the control variables chosen for this study were: 

• jet exhaust area (u^) 

• fan inlet guide vanes (ug) 

• rear compressor variable vanes (u^) 

• main burner fuel flow (u^) 

The engine state variables chosen were: 

• fan turbine inlet temperature (x^) 

• main burner pressure (xg) 

• fon speed (x-) 

• high compressor speed (x^) 

• afterburner pressure (x^) 


b 
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These independent variables are sufficient to establish the F401 operating condition 
and to define the dynamic path of the engine. The engine output variables were 
selected after consideration of engine steady-state operating requirements. For 
efficient steady-state performance the positions of the fan inlet guide vanes and 
the rear compressor variable vanes are scheduled as functions of fan speed, N^, 
and high compressor speed, Ng, respectively. It is therefore desirable that, in 
the steady state, the vane positions be at their scheduled values. This was accom¬ 
plished by (l) defining the fan inlet guide vanes and the rear compressor variable 
vanes as the first two output variables, yi and y2, respectively, and (2) commanding 
that yp and y2 achieve their steady-state scheduled positions as a function of 
steady-state Np, 1respectively. This procedure allowed transient vane operation 
to be optimised wnile insuring proper steady-state positioning. The remaining 
engine output variables considered were: 

• thrust (y^) 

• high turbine inlet temperature (y^) 

• fan corrected airflow (y^) 

• fan stability margin (y^) 

• compressor stability margin (y^) 

Table I lists values of the steady-state engine variables for sea-level static 
conditions at five power-lever angle (PLA) design points: PLA = 20, 35, 47, 60 and 
73 deg. For convenience, all engine parameters except fan and compressor vanes and 
fan and compressor stability margins have been normalized to 1.0 at FLA. = 73 deg. 

The vane positions are defined as ration of their maximum positions, and the two 
stability margins are given as ratios of one (a smaller value indicates reduced 
stability margin). These data are also shown in Fig. 1, where linear interpola¬ 
tion has been used between design points to define the steady-state values as 
functions of power-lever angle. The 20-deg PLA is the engine idle condition while 
the 73-deg PLA represents the military rating condition (maximum nonafterburning 
thrust). 


The next step in the synthesis procedure is to linearize the engine dynamics 
at a set of points along the steady-state operating line. Linear dynamics at a 
steady-state operating point can be represented by 


6 x = A6x + B6u 
6y = C6x + D6u 


( 1 ) 
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where the vector Sx represents perturbations in the five engine state variables, 
the vector Sy represents perturbations in the seven output variables, and the 
vector 8xi represents perturbations in the four control variables. The constant 
matrices A, B, C and D consist of appropriate partial derivatives of the engine 
dynamic response evaluated at the given operating point. Linearized engine dynam¬ 
ics at PLA = 20, 35, 47, 60 and 73 deg, calculated in a previous program (Ref. 1), 
were used here for control synthesis. The calculated values of A and B at each of 
the five control design points are presented in Table II where the control, state, 
and output variable ordering is as previously defined. The corresponding C and D 
matrices are presented in Table HI. The matrix elements were determined by apply¬ 
ing a system identification technique to input-output-state data generated by a 
nonlinear dynamic computer simulation of the F401 engine (Ref. 3). Details of the 
system identification procedure used are contained in Ref. 1. 

Four integrators, one for each control input, were also added to the linear 
engine model of Eq. (1). These integrators are used in the optimal controller to 
insure zero steady-state errors in commanded engine response. Augmenting the 
linear engine dynamics with the integrators results in the following equations 
for the controlled plant: 

6x = A6x + B6u 

6u = 6w (2) 

6y = C6x + D6u 

The linear engine model with the integrators is shown in Fig. 2(a). 


Calculation of Optimal Feedback Gains 

Optima], linear feedback gains were computed at the five selected control design 
points by minimizing the following performance index: 


» r 

= ij I q 1 (lGV) 2 + ^(RCW) 2 + q 3 (F) 2 + q i+ (T IT ) 2 + q^(FSM) 2 + (^(CSM) 2 


+ r l(A j ) + r 2 (lGV) + r 3 (RCW) 2 + r 4 (w f )‘ 


dt 


where q^, r^ are weighting factors and 

IGV = perturbed fan inlet guide vanes 

RCW = perturbed rear compressor variable vanes 


(3) 
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F 5 

perturbed 

*?_, it 

i I? - 

perturbed 

FSM = 

perturbed 

CSM i 

perturbed 

A i - 

perturbed 

w f k 

perturbed 


thrust 

high birbine inlet temperature 
fan stability margin 
compressor stability margin 
jet exhaust area 
main burner fuel flow 


The act notation denotes time rate-cf-change of a parameter. The solution tc the 
linear optimal control problem, i.e., the mi nimi zation of J (Ref. 4), leads to the 
standard optimal control structure shown in rig. 2(b). The constant matrices 3^ 
and Gp (Fig. 2(b)) are the optimal closed-loop feedback gain matrices. The gain 
matrix M has been added to the controller tc permit the consideration of command 
inputs 6z. 


Tc insure zero steady-state errors between desired and actual output, an 
integral control structure (Fig. 2(c)) was derived from the standard form. A 
method for converting from the common linear optimal structure to an integral 
control which retains the desired optimal control characteristics was developed 
in Ref. 1. The integral control gain matrices H and L are determined from the 
system matrices A, B, C and D and the standard optimal gain matrices G^ and Gg. 
Since there are four independent engine control variables, the steady-state values 
of four engine output variables can be specified independently. The vector v 
(Fig. 2(c)) represents those four engine output variables whose steady-state values 
are to be specified. 


Nonlinear Feedback Control Synthesis 


After the H and L matrices (Fig. 2(c)) were computed at each control ierirn 
point, linear interpolation was used to define the matrices as functions of nigh 
compressor speed, N2, between design points, i.e., H = H^), L = L(K^). If H 
and L are defined in this manner, the optimal control is 


u*(t) = 


fx(t) 0 

J L(N 2 ) dx + j H(Np ) fv(r) -z(t) 1 dT + u*(0) 
x(0) QJ 


(4) 
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The vectors v and z are 


= fan inlet guide vanes 
v 2 = rear compressor variable vanes 
v^ = thrust 

v^ = high turbine 3nlet temperature 


Z 1 

z 2 

s 3 

z 4 


commanded v 7 
commanded v 0 
commanded v^ 
commanded v^ 


The functions z^, Zg, z 5 , z^ were coordinated Such that the engine was commanded 
from initial to final conditions both of which were always on its steady-state 
operating line (as defined in Fig. 1) when operating point was changed. The 
resulting nonlinear closed-loop control structure is shown in Fig. 3. For con¬ 
venience, the engineering symbols used in Fig. 3 fur the engine variables denote 
actual rather than perturbational values. 


Characteristics of the Optimal hfciltivariable 
Controller Developed for the F401 Engine 

The performance index weightings a^, r^ (see Eq. (3) for definitions) were 
selected at each control design point so that the resulting controller feedback 
gains (1) produced a rapid engine acceleration and (2) avoided exceeding maximum 
allowable control actuation rates, maximum control positions, and minimum allow¬ 
able fan and compressor stability margins (stability must be greater than zero). 
Table IV presents maximum and minimum limits on F401 control positions and rates. 
Table V lists a set of weightings for which engine response with the optimal con¬ 
troller will be shown. The higher thrust weightings, q^, at PLA = 47, 60 and 
73 deg (Table V) result in a rapid thrust response during that part of the tran¬ 
sient where compressor and fan stability margins pose no problem. The large weight¬ 
ings on high turbine inlet temperature, q^, at PLA = 20 and 35 deg result in a very 
rapid temperature response, which in turn helps to provide a rapid engine accelera¬ 
tion. The even larger weighting on temperature at PLA = 47 deg regulates the 
temperature response when temperature is close to its maximum commanded value, 
and thereby prevents overtemperature. 

The weightings on rate-of-change of fan stability margin, q^, at PLA = 47 and 
60 deg (Table V) permit rapid accelerations at high power levels without surging the 
fan (e.g., engine accelerations from 80 to 100 percent thrust). Such a weighting 
is not necessary at the idle condition (PLA = 20 deg) because the large weighting 
there on rate-of-change of jet exhaust area, r-^, protects against fan surge. Also, 
it was found that if r^ were chosen too small at PLA = 20 deg, the jet area would 
respond to a large step in PLA by immediately decreasing to its bottom position 
limit. Jet area would then remain there during most of the resulting transient 
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before increasing to its steady-st *te value. This presented a problem, since the 
movements of the remaining three controls were being defined on the basis of a jet 
area which was free to move in either a positive or negative direction. Consequent¬ 
ly, the movements of these three controls (fuel flow, fan inlet guide vanes, rear 
compressor variable vanes) became uncoordinated with the jet area me/ament, end poor 
closed-loop performance resulted. The larger weighting on rate-of-change of jet 
area at PLA. = 20 deg thus serves the dual purposes of (1) maintaining coordination 
between all four control variables by preventing the jet area from operating on its 
bottom limit for a long period of time, and (2) preventing fan surge by not allow¬ 
ing the jet area to immediately decrease rapidly after a step PLA. change. 

The high weighting or. rate-of-change of compressor stability margin, qg, at 
PIA = 20 deg (Table V) insures adequate compressor margin at the beginning of the 
transient, where compressor stall is most likely to occur. Increased levels of 
compressor margin are also obtained b^ not allowing the rear compressor variable 
vanes to move rapidly away from their steady-state schedule during the initial 
stages of the transient. This is the purpose of the high weightings at PIA = 20 
and 35 deg on rate-of-change of compressor vanes, r^. 

Feedback gain matrices H and L which result for the. weightings of Table V are 
presented in Table VI. Recall that intermediate values of H and L were computed 
for the optimal controller by interpolating linearly with r^. Consequently, H 
and L can he regarded as continuous functions of No between the end points shown 
in Table VI. 


Comparison of F401 Engine Response with Opti m al 
Multivariable and Conventional Controllers 

The optimal multivariable feedback controller (Fig. 3) having the parameters 
described in the previous section was implemented on the UARL nonlinear computer 
simulation of the F401 engine (Ref. 3)» Engine response for an acceleration command 
from idle to military thrust, i.e., a step change in PLA. from 20 deg (9 percent 
thrust) to 73 deg (100 percent thrust), was calculated. The results are shown in 
Figs. 4 through 10. Figures 4 through 8 show time responses for the following 
engine variables, respectively: thrust, high turbine inlet temperature, compressor 
and fan stability margins, fan and compressor speeds, airflow, exhaust area and fuel 
flow. Changes in the positions of the fan inlet guide vanes and rear compressor 
variable vanes are plotted versus fan and compressor speeds in Figs. 9 and 10, 
respectively. Fan and compressor variable geometries are presented as functions 
of spool speeds to more clearly depict the relationship betwton transient and 
steady-state vane positions. Comparative engine response data for an F401 con¬ 
ventional controller, which was designed to provide rapid engine acceleration, 
are also shown in Figs. 4 through 10. Conventional response data were generated 
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by tne Florida Research and Development Center of Pratt & Whitney Aircraft using 
an FUOl engine deck similar t' that provided for UARL. However ; there are slight 
differences in the decks which cause the small mismatch evident between the two 
sets of results at the steady-state operating points. 

Thrust response for the optimal controller is approximately 30 percent more 
rapid than for the conventional controller (Fig. h). That is, the optimal cover oiler 
thrust response reached the 93-percent pcint in 2.32 sec, while conventional thrust 
response reached the same point in 3*^0 sec. The high turbine inlet temperature 
response for the optimal control was also significantly faster than the conventional 
temperature response (Fig. 5(a)). This difference contributed to the thrust improve¬ 
ment fer the optimal system. Minimum levels of compressor and fan Suability margins 
were similar for the two types of control (Figs. 5(b) and 5(c)), as were tic fan and 
compressor speed and the airflow responses (Fig. 6). 

Comparing t ime histories of the four engine control variables for the different 
controllers following the step PLA. command (Figs. 7 through 10) provides insight 
into how the optimal controller achieved its rapid acceleration. The optimal con¬ 
troller moved jet exhaust area (Fig. 7) and fan inlet guide vanes (Fig. 9) in a 
manner distinctly different from the corresponding conventional control movements. 

For the optimal controller, jet exhaust area (Fig. 7) did not move far from its 
initial position until the fan stability margin (Fig. 5(c)) had begun to increase 
(about 1.1 sec into the transient). At that time, jet exhaust area decreased rapid¬ 
ly to its bottom position limit, and this was followed by a rapid increase to its 
final steady-state value. With the conventional controller, exhaust area increased 
rapidly, remained constant, and then decreased to its steady-state value (Fig. 7). 

Fan inlet guide vanes, under optimal control, led the programmed steady-state 
schedule*, but they lagged the steady-state schedule with the conventional controller 
(Fig. 9)* These differences between the way in which optimal and conventional con¬ 
trollers change the F401 variable engine geometries are similar to those previously 
reported for high-thrust (80 to 100 percent) engine accelerations (fief. 2). 

Finally, engine fuel flow input for the optimal controller was much more rapid 
than for the conventional controller (Fig. 8). This difference in engine fuel flow 
resulted in the .core rapid temperature response of the optimal system which was 
described previously (Fig. 5(a))* The optimal controller was able to rapidly input 
fuel flow, but, by the proper coordination of all engine controls, avoid overtem¬ 
perature and stability problems. 


*The steady-state schedule adjusts fan and compressor vane position with fan and 
compressor speeds, respectively, to insure efficient steady-state engine operation. 
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SYNTHESIS CF AN OPTIMAL INTEGRATED 
INLET-ENGINE FEEDBACK CONTROL SYSTEM 


In this section the design of an integrated feedback controller for the FlOl 
turbofan engine model and an internal comrression supersonic inlet model is de¬ 
scribed. The design objectives end conditions are described initially. Inlet- 
engine input, state and output variables appropriate for synthesizing feedback cen¬ 
tre-. are defined next, and the coupling of inlet-engine dynamics is discussed. The 
linear optical control methods used tc compute feedback giins are then outlined, 
and, finally, propulcion syszem response to simulated afterburner ignition for 
integrated and separate inlet and engine controllers is compared. 


Design Objectives and Conditions 

The objective cf this effort was tc design an integrated controller which uses 
cross-coupling feedback paths between the inlet and the engine to provide improved 
inlet-engine dynamic performance, linear optimization method'- were used t r compute 
integrated control gains, and also to compute the inlet control and engine control 
gains separately. The latter results were used tc provide comparative inlet-engine 
response for non-integrated controls. Design criteria upon which control synthesis 
was based were (.1; regulating the inlet normal shoes position, (2) regulating the 
throat Much number, and (3) maintaining smooth far. airflow. In addition, the con- 
trcl designs were to allow only small variations in engine stabilit” margins and 
high turbine inlet temperature sc as to insure safe closed-loop system response. 

The effectiveness cf the controllers in satisfying these design criteria was 
evaluated for simulated engine afterburner ignition. 


Supersonic Inlet Model 

An internal compression supersonic inlet model was selected for this study 
because this type cf inlet provides the greatest potential, for improved dynamic 
performance. The internal compression inlet operates at its greatest efficiency 
(highest pressure recovery) when the threat Mach number approaches one, and the 
normal shock is located very near the inlet throat. However, the no rm a l shock 
will be expelled from the inlet if the normal shock reaches the throat. This 
phenomenon is termed "inlet unstart" and results in very inefficient pressure 
recovery as well as distorted airflow to the engine. An integrated inlet-engine 
controller may result in significantly increased pressure recovery by permitting 
safe operation closer to the unstart condition. 
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A mathematical model of an internal compression supersonic inl et linearized 
about a flight condition of 40,000 ft, Mach 2.2 was constructed by the Hamilton 
Standard Division of United Aircraft Corporation. This flight condition was chosen 
because it represents a military-specified operating point at which aircraft/engine 
performance is evaluated. Inlet compatibility with the F401 engine was achieved b> 
designing the inlet for ?401 fan airflow and fan airflow Mach number at the selected 
flight condition. The linearized inlet model is shewn in Fig. 11. The inlet con¬ 
trols selected were throat area (u^p/j bypass area (u^?), and fan corrected airflow 
(uq). Fan corrected airflow (an output of the engine) is the coupling variable 
from the engine to the inlet. The outputs of the two integrators in the inlet model 
(Fig. 11) were defined to be the inlet states, xjp and Xj 2 * The outputs chosen for 
the inlet were throat Mach number (yip), normal shock position (yjp)> and fan inlet 
total pressure (yj^). 

The A, B, C and D constant matrices of Eq. (1) for the inlet model ar 2 given 
in Table VII. Table VIII lists values of the steady-state inlet and engir.e param¬ 
eters at the 40,000 ft, Mach 2.2 design point. For convenience, threat area, fan 
inlet total pressure, and afterburner fuel flew were normal,red to 1.0 at the 
design flight condition. Also, normalized values cf bypass area and normal shock 
position were defined to be 0.11 and 0.16, respectively, at the design flight condi¬ 
tion. These values for bypass area and normal shock position correspond to their 
actual values when they are expressed in ft 2 and ft, respectively. Note that a 
variation cf -0.l6 in the normal shock position represents the iilet unstart condi¬ 
tion. Normalization of the other engine parameters was identical to that defined 
previously in the discussion of the engine controller sy-thesis. 


Identification of Coupled Inlet-Engine Dynamics 

Linearized engine dynamics at the 40,000 ft, Mach 2.2 flight condition were 
identified first as a prelude to interfacing the inlet model of Fig. 11 with tne 
F401 engine model. Engine state variables were as previously defined; however, 
the engine control vector was augmented by afterburner fuel flow and fan inlet 
total pressure. Note that perturbed fan inlet total pressure is also an output 
of the inlet (_,ee Fig. xl); hence, it represents the coupling variable from the 
inlet to engine. The engine output variables considered here were thrust, high 
turbine inlet temperature, fan corrected airflow, fa., stability margin, and com¬ 
pressor stability margin. The A, B, C and D constant matrices for this operating 
point were determined by applying a system Identification technique to input-output- 
state data generated by the F401 computer simulation. Details of the system identi¬ 
fication procedure used are contained in Ref. 1. The computed values of the A, B, 

C and D engine matrices with input, state, and output ordering as defined in this 
section ere presented in Table IX. 
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The next ste- in the identification of coupled inlet-engine dynamics was to 
represent the ccznbined inlet-engine system in the form of 33a. (1). The overall 
inlet-engine model was formed by combining the previously developed inlet and 
engine models. The combined inlet-engine states and outputs were the previously 
defined engine and inlet states and outputs, respectively. Controls for the inlet- 
engine system consisted of the previously defined inlet and engine controls, with 
the exception of the coupling variables. Coupling variables represent parameter., 
which are internal to the overc system. The coupling variable from the inlet to 
the engine was fan inlet total pressure and the coupling variable from the engine 
to the inlet was fan corrected airflow. The resulting combined system A, 3, C 
and D matrices as well sis control, output, and state variable ordering are presented 
in Table X. 


Synthesis of Integrated Inlet-Engine Control 
Optimal linear feedback gains were computed by minimizing the performance inde: 


j = 


cc 

£ / q^Trj ,) 2 + q 2 (w aix ) 2 + Q /ISM ) 2 + q^CSM ) 2 + q^Nj ,) 2 + q ^) 2 


+ ayO-fo) + qe(x s ) ■ r^f + r 2 (lGV) + ^(RCW) + r u (w f ) 
+ r 5 ^ab)^ + r 6 ( V" + r 7 ( *BY r 


where q^, 


r^ are 


w A 
anr ~ 


Np 


N r 


4 
4 

Mjh = 
x 4 
4 
4 

^BY = 


w 


ab 

Ap H 


weighting ' ^rs, and 
perturbe' fan corrected airflow 
perturbed fan speed 
pertrrbed high compressor speed 
perturbed throat Mach number 
perturbed normal shock position 
perturbed afterburner fuel flow 
perturbed throat area 
perturbed bypass area 


(5) 
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while the other symbols are as previously defined. The performance index includes 
all of the input rates to account for and satisfy actuation rate limits. The remain- 
in*, parameters are included because their response is critical in evaluating closed- 
loop system perfor man ce. The structure of the resulting multivariable controller 
is shewn in Fig. 2(b}. Because only the regulator problem was considered (6z = 0 
in Fig. 2(b)), the standard structure of Fig. 2(b) was not transformed to the 
integral structure of Fig. 2(c). 

Before the weighting factors were chosen, acceptable propulsion system response 
to a simulated afterburner ignition was defined. Hie design —it aria for inlet per¬ 
formance specified that variations in threat Mach number and normal shock position 
be minimized. For satisfactory engine performance the variations in fan corrected 
airflow, fan stability margin, compressor stability mar 6 ln and high turbine inlet 
temperature were also to be minimized. In addition, the engine was to be maintained 
near its steady-state operating point since this point represents efficient engine 
operation for the given flight condition. Ibis was accomplished by permitting only 
small transient variations in fan and high compressor speeds. 

The performance index weightings were selected so that the resulting controller 
feedback gains produced satisfactory inlet-engine response tc a disturbance equiva¬ 
lent to that from an afterburner light. This disturbance consisted cf a step change 
in afterburner fuel flaw. A very large weighting factor, r^, was used to insure 
that afterburner fuel flow had a negligibly small rate-of-decrease to zero from 
its stepped value. Table XI lists the set of weightings and Table XII the con¬ 
troller feedback gains for which inlet-engine response 'will be shown. The inte¬ 
grated controller is comprised of engine feedback gains, inlet feedback gains, and 
cross-coupling gains between the inlet and engine. The gain matrices (G^ and G2, 

Fig. 2(b)), have therefore been partitioned as 


G 1 ~ 

>_ 

w 

G 2 = 

" G ?v ' 
| 

G 25I 


. G 1IE 

W 

'S* 


_ G 2IE 1 

G 2I 


The 5x5 Gj_g and G2E matrices represent feedback paths to engine control rates 
frc.a engine states and control inputs, respectively. The feedback gains to inlet 
control rates from inlet states and controls are given by the 2x2 gain matrices 
G11 and G2j } respectively. The 5x2 gain matrices G^ei and G2 ej re P reyen t the 
coupling gains from the inlet states and controls, respectively, to the engine 
control rates. The coupling gains from the engine states and controls to the .inlet 
control rates are the 2x5 gain matrices G^-j-g and G2135, respectively. 


14 



N9H620-2 


Synthesis of Separate Inlet and Engine Controllers 

Separate inlet and engine controllers were also designed using modern control 
theory. The separate inlet controller used only inlet information to achieve satis¬ 
factory inlet performance, and the separate engine control u^ed only engine informa- 
tion. Designing the engine ana inlet controls separately ignored the known cross¬ 
coupling between engine and inlet variables. However, there was no specification 
that a particular control variable was to control a particular output variable. 
Cross-coupling paths within the engine alone and within the inlet alone were used. 

Feedback gains for the inlet control were computed by minimizing 


ac r 

J = t f ^(Mpg) + Qe(x s ) * r £^H^ + r 7^Y^ + r 8^ w air^ 

0 1 


dt 


(?) 


and the gains for the engine control were computed by minimizing 

J = if [ q 1 (^ir) 2 + a -£( w air^ 2 + a -3^ FSM ) 2 + <V CSM ) 2 + < 3 5 ( h f) 2 + 


( 8 ) 


+ r-jU^) 2 + r 2 (lG7) 2 + r^RCW) 2 + r 4 (w f ) 2 + r^w^)' 


dt 


where the symbols are as previously defined. Recall that for the inlet model alone, 
perturbed fan corrected airflow was considered an input (see Fig. 11), The control 
optimization procedure therefore theoretically called for direct fan airflow control 
to regulate the inlet. However, this is inconsistent with propulsion system design. 
That is, fan airflow can only be controlled indirectly, e.g., by changing system 
geometries or fuel flows. To reconcile the model mathematics with the actual pro¬ 
pulsion system, a very large weighting, on perturbed airflow was used. The 
large weighting resulted in negligibly small feedback gains to airflow and thereby 
avoided the implications of direct airflow control. 


Weighting factors were chosen for the performance indices of Eqs. (7) and (8) 
which were identical to the weighting factors in the integrated control performance 
index of Eq. (5). The gain matrices Gp and Gp for the inlet and engine models, 
when the inlet controller and engine controller were designed separately, are given 
in Table XIII. Note that the gain matrices resulting from separate engine and inlet 
control synthesis are approximately the same as the integrated control gain matrices 
(Table XII), except that the inlet-engine cross-coupling feedbacks are zero for 
separate control synthesis. 
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Coaroartson of Inlet-Engine Response 
with Separate and Integrated Controllers 

Inlet-engine node! response with the integrated controller was calculated for 
a simulated afterburner ignition consisting of a 10-percent step in afterburner fuel 
flow. The effects of this step are considered to be representative of the distur¬ 
bances caused by afterburner lightoff. The closed-loop time responses for the 
integrated controller are shown in F_g. 12. For the given step in afterburner fuel 
flow, a C percent increase in steady-state thrust resulted with less than 1-percent 
maximum variation in fan stability margin and less than 0.5-percent maximum varia¬ 
tion in far* corrected airflow, compressor stability margin end high turbine inlet 
temperature. The changes in the steady-state values of the fan and high compressor 
speeds caused by simulated afterburner ignition were negligible. The maximum varia¬ 
tion in throat tech number was less than 0.002 and the maximum variation in normal 
shock position was 0.0053* 

Comparative closed-loop time responses are also bhcwn in Fig. 12 for the con¬ 
troller which had been designed by considering the engine and inlet separately. 

The engine performance remained essentially unchanged from that which resulted for 
integrated control. There was a minor improvement in throat tech number variation; 
however, the maximum variation in normal shock position more than tripled from 
-0.0017 for integrated control to -0.0053 for separate control (recall that a varia¬ 
tion of -0.16 would cause inlet unstart). This tripling of the normal shock posi¬ 
tion variation represents a significant deterioi ation in dynamic performance due to 
separate inlet and engine control. These results also imply that for integrated 
control it may be possible to operate closer to the inlet unstart condition and 
thereby achieve greater inlet pressure recovery. 

As noted before, the gems of the integrated controller are very similar to 
the gains of the separate inlet and engine controllers, except that the integrated 
controller contains inlet-engine cross-coupling terms. To obtain further insight 
into the influence of these cross-coupling terms, the inlet-engine model response 
was calculated for the integrated controller with cross-coupling gains zeroed. 

Results showed that inlet-engine performance with the cross-coupling terms removed 
from the integrated controller was almost identical to system performance with the 
separately designed controller. This indicates that the improved xnlet performance 
with the integrated control was a direct result of the cross-coupling feedback gains. 

To determine whether a controller without inlet-engine cross-coupling gains 
could he designed to perform as well as the integrated controller, the normal shock 
position weighting, qg, in the inlet control performance index (Eq. (7)) was in¬ 
creased. Recall that it was the normal shock whose response deteriorated the most 
with the separate controllers. As the change in inlet normal shock position was 
penalized more heavily, the observed variation in normal, shock position for the 
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separate controllers decreased, but the variation in throat Mach number increased. 
The engine response remained essentially unchanged. When the weighting (qg) was 
increased to the point where the throat Mach number variation was the same as that 
for the integrated control (recall that it had been less for the separate controls 
with the original ag weighting), the normal shock position variation was approxi¬ 
mately double that of the integrated control. The weighting was then further 
increased until the variation in normal shock position was essentially the same as 
that for the integrated control. Inlet response for this controller is compared 
with inlet response for the previously developed integrated controller in Fig. lj. 
The throat Mach number variation is six times larger with the separately designed 
controller than with the integrated controller. The bypass area excursion also 
increased by a factor of four while the throat area excursion increased by a factor 
of six. Both inlet actuators were required to move twice as fast with the separate 
controllers as with the integrated controller. Because the engine control remains 
unchanged, and the computed inlet variations due to the afterburner ignition have 
little effect on the engine, the engine response was identical to that presented 
in Fig. 12. The inlet weighting factors and inlet gain matrices (G^j and G21) 
for this high qg weighting controller (03 = 50,000) are shown in Table XIV. The 
t.jgine control was not changed so that its parameters remain as given in Table XIII. 

The results discussed here demonstrate that inlet response to an afterburner 
ignition can be improved by an integrated controller. The optimal control theory 
identified fan corrected airflow as a coupling variable and defined cross-coupling 
gains to provide improved inlet regulation of the effects of a downstream distur¬ 
bance. The cross-coupling feedbacks from engine to inlet modulated the inlet con¬ 
trols so as to maintain normal shock position with a smaller variation in throat 
Mach number end less control effort than the separately designed controller. The 
engine response to inlet distortion was not investigated. However, the computed 
coupling gains from the inlet to engine verify that the fan inlet total pressure 
coupling variable was identified by the control synthesis procedure. These gains 
could be used to modulate the engine inputs and improve engine performance to an 
upstream disturbance. 
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CURRENT UARL FRCX®AM 


The current UARL program sponsored by ONR is focused on the development of 
parameter identification and adaptive control for advanced technology propulsion 
systems*. Identification of engine parameters is the key to applying the control 
optimization method described in this report to a real rather than simulated engine. 
That is, the piecewise-linear/piecewise-optimal approach tc multivariable control 
synthesis requires that the engine stability and control derivatives be known (A, 

B, C and D matrices). Accurate determination of turbofan engine stability and 
control parameters is difficult, because (1) not all engine variables can be sensed, 
(2) data from sensed engine variables are noise-corrupted, and (3) the parameter 
values change with altitude, Mach number, and power lever angle. Algorithms which 
account for these difficulties and provide valid estimates of the engine dynamics 
are currently not available. 

Modern (optimal) filtering methods are applicable to realistic jet engine 
estimation problems, e.g., operation on noise-corrupted engine input-output data 
and incomplete sensing of engine variables. Optimal filtering theory is being 
employed in the current UARL studies to estimate those parameters which are criti¬ 
cal in assessing engine performance. Also, on-line computation of engine parameters 
provides the data necessary to develop and implement a self-adaptive control algo¬ 
rithm which would optimize propulsion system performance in flight. The feasibility 
of using these newly estimated parameters as the basis for on-line adaptive control 
optimization will be evaluated in the UARL program. An algorithm for adaptive 
engine performance control will be developed using thrust specific fuel consumption 
as the performance criterion for the optimization process. 


*In a related program being funded by P&WA, UARL and P&WA engineers are evaluating 
hardware and software requirements involved in implementing the optimal multivariable 
controllers developed under ONR sponsorship. A PDP-11 digital computer is being 
coded to simulate the optimal multivariable controller ir. real time. An F401 engine 
model, defined by the parameters of Tables II and III, is being programmed on a 
Beckman analog computer. The resulting real-time hybrid simulation will be used 
(1) to study system accuracy and closed-loop stability, (2) to investigate digital 
control cycle time requirements, and (3) to resolve interface problems associated 
with a discrete controller and a continuous F401 engine. 
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A 


A BY 


%’H 

B 


C 

CSM 

D 

deg 


E 


F 

FSM 

G 


G 1 

G 1E 

g iei 

G 1I 

^IE 

G 2 

G 2E 

G 2EI 


LIST OF SYMBOLS 

Constant n x n matrix used to describe linearized system dynamics 
Perturbed bypass area, normalized 
Perturbed jet exhaust area, normalized 
Perturbed throat area, normalized 

Constant a x m matrix used to describe linearized system dynamics 

Constant p x n matrix used to describe linearized system dynamics 

Perturbed comprensor stability margin, normalized 

Constant p x m matrix used to describe linearized system dynamics 

Degrees 

Constant k x 7 matrix relating system output tc commanded output 

Perturbed thrust, normalized 

Perturbed fan stability margin, normalized 

Standard optimal closed-loop feedback gain matrix 

m x n partition of the matrix G 

5x5 partition of the inlet-engine gain matrix G]_ 

5x2 partition of the inlet-engine gain matrix G^ 

2x2 partition of the inlet-engine gain matrix G-^ 

2x5 partition of the inlet-engine gain matrix G]_ 
m x m partition of the matrix G 
5x5 partition of the inlet- ,ngine gain matrix G 2 
5x2 partition of the inlet-engine gain matrix G 2 
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r* 

U2i 

•hiE 

H 

IGF/ 

i 

•T 

L 

M 

E 

N C 

N F 

K 1 

r 2 

n 

PLA 

p 

% 

RCW 

r i 

t it 


LIST 0? SYMBOLS (Continued) 

2>:2 partition of the inlet-engine gain matrix Gg 
2x5 partition of the inlet-engine gain matrix Gg 
4 x 4 optimal integral feedback gain matrix 
Perturbed fan inlet guide vanes, normalized 
General subscript 
Performance index 

4x5 optimal integral feedback gain matrix 

m x m feedforward gain matrix 

Dimension of system control vector u 

Perturbed throat Mach number 

Ferturbed high compressor speed, normalized 

Perturbed fan speed, normalized 

Fan speed, normalized 

High compressor speed, normalized 

Dimension of system state vector x 

Power-lever angle 

Dimension of system output vector y 
Performance index weighting factor — a scalar 
Perturbed rear compressor variable vanes, normalized 
Performance index weighting factor — a scalar 
Perturbed high turbine inlet temperature, normalized 
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LIST OF SYMBOLS (Continued) 


t Time, sec 

u in x 1 control vector 

u* m x 1 optimal control vector 

i^ 1 component of u 

Ujj inlet control variable 

v 4x1 vector consisting of those elements of the output vector y 

whose steady-state values it is desired to specify 

v^ i^' n component of v 

w m x 1 control rate vector — w = du/dt 

w ab Perturbed afterburner fuel flow, normalized 

"air Perturbed fan airflow, normalized 

w f Perturbed main burner fuel flow, normalized 

x n x 1 state vector 

x^ i^h component of x 

x s • Perturbed normal shock position, normalized 

Xjj inlet state variable 

y p x 1 output vector 

y^ i^ h component of y 

y-j-j j a inlet output variable 

z m x 1 command inpv>t vector 

i^* 1 component of z 
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LIST OF SYMBOLS (Concluded) 

t Variable of integration 

(*) Time derivative of the quantity in parentheses 

d( ) Differential of the quantity in parentheses 

6( ) Perturbational value of the quantity in parentheses 

A Equals by definition 
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TABLE I 

HC8HMJZBD STEADY-STATE ENGINE PARAMETER 
AS A FUNCTION CF POWER LEVER ANGIE 


Engine Parameters 

Power Lever Angle, 

PLA-deg 

Type 

Parameter 

20 

35 

47 

bC 

73 


Thrust 

0.09 

0.35 

0.52 

0.72 

1.0 


High Turbine Inlet Temperature 

0.56 

0.73 

0.81 

0.90 

1.0 

Output 

Fan Airflow 

0.3b 

0.62 

0.73 

0.85 

1.0 


Fan Stability Fargin 

0.10 

0.14 

0.18 

0.12 

0.12 


Compressor Stability Kargin 

0.16 

0.18 

0.18 

0.19 

0.19 


Fan Turbine Inlet Temperature 

0.57 

0.72 

CO 

• 

0 

0.90 

B9 


Kain Burner Pressure 

0.21 

0.45 

0.59 

0.77 


State 

Fan Speed 

0.49 

0.74 

O.83 

0.91 

1 


High Compressor Speed 

0.69 

O.83 

0.88 

0.93 

1.0 


Afterburner Pressure 

0.40 

0.55 

0.66 

0.82 

1.0 


Jet Exhaust Area 

O.98 

O.98 

0.98 


181 

Control 

Fan Inlet Guide Vanes 

-0.50 

-0.50 

-0.50 


EE9 

Rear Compressor Variable Vanes 

-1.11 

-0.39 

-0.17 


mm 


Main Burner Fuel Flow 

0.12 

0.33 

0.48 


1.0 
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TABLE ZL 

IDENTIFIED VALUES G? LINEARIZED ENGINE DYNAMICS - 
A AND 3 MATRICES AS A FUSCTIC8? CF POWER LEVER ANGIE 


PLA-deg 

Matrix 

Matrix Elements 



-56.450 

19.387 

2.403 

-48.947 

17.867 



8.066 

-71.982 

2.184 

47.511 

-1.2V* 


A 

0.123 

4.135 

-1.672 

-0.50G 

-2.463 



0.222 

3.912 

-0.098 

-2.791 

-0.514 

20 


-0.877 

5-472 

1.309 

-2.938 

-9.303 


0.626 

0.005 

-1.24? 

103.080 





-0.183 

0.001 

1.362 

3.300 



B 

-0.011 

-0.085 

-0.011 

0.030 




0,004 

0.003 

-0.097 

-0.131 




-0.917 

0.017 

-0.147 

0.198 




-64.848 

12.421 

-I5.U.3 

-36.828 

2.202 



23.201 

-68.8^8 

24.036 

66.740 

7.405 


A 

0.448 

4.983 

-2.911 

-0.766 

-2.742 



1.253 

2.071 

-0.473 

-2.881 

-0.777 

35 


-1.720 

7.956 

-0.538 

-3.743 

-7.617 


0.808 

-0.144 

-4.129 

56.602 





0.020 

0.260 

5.866 

-2.124 



B 

-0.063 

-0.132 

-0.055 

-0.316 




-0.074 

-0.011 

-0.180 

-0.733 


• 


-2.856 

-0.258 

-0.589 

-0.397 




-57.096 

3.613 

-10.211 

-5.481 

-2.715 



19.832 

-72.340 

30.295 

40.972 

15.327 


A 

0.660 

4.496 

-3.601 

-0.011 

-2.808 



1.326 

2.313 

-0.809 

-3.032 

-0.821 

1 ~ 


0.882 

0.703 

2.922 

1.471 

-4.596 

*+( 


1.017 

-0.553 

-3.941 

39-792 




-0.125 

l.4l6 

7.888 

4.181 



B 

-0.077 

-0.316 

-0.031 

-0.382 




-0.088 

-0.033 

-0.253 

-0.565 




-3.563 

-0.149 

-O.097 

-0.785 



(Continued) 
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II (Concluded) 


Matrix Elements 


-21.855 

4.762 

8.122 

4.711 

-31.287 

9.385 

15.460 

4.615 

4.729 

-3.880 

-0.156 

-3.095 

1.991 

-0.828 

-2.524 

-0.817 

-2.882 

4.688 

3.235 

-3.436 

0.370 

-0.273 

27.141 


0.598 

4.171 

9.866 


-0.394 

-0.036 

-0.362 


-0.065 

-0.367 

-0.513 


-0.069 

0.339 

0.396 


-9.303 

12.037 

-2.398 

-1.254 

-38.762 

-4.221 

28.480 

14.729 

2.287 


-1.546 


1.062 

-0.729 

-2.150 

-0.624 

-8.814 

-0.167 

7.477 

1.099 

0.546 

-0.813 

17.095 


1.341 

7.737 

8.641 


-1.176 

-0.4i6 

2.034 


-0.024 

-0.55" 

-0.378 


0.874 

1.617 

0.223 
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(TABLE m 

IDEI<TI?ISD VAUJES 0? LINEARIZED ENGINE DYNAMICS - 
C AND D MATRICES AS A FUNCTION OF POWER IEVER ANGLE 






PIA-deg 

2-Jatrix 

Matrix Elements 



0 

0 

0 

0 

0 



0 

0 

0 

0 

0 



- 0.006 

-0.025 

- 0.013 

0.002 

2.118 


C 

1.057 

0.471 

- 0.028 

- 0.028 

- 0.228 



- 0.023 

0.075 

1.103 

- 0.030 

-0.304 



-0.053 

0.314 

0.329 

- 0.109 

-1.546 

20 


-0.244 

-4.114 

- 0.026 

4.415 

- 0.023 


0 

1.0 

0 

0 





0 

0 

1.0 

0 




0.095 

-0.002 

0 

0.009 



D 

-0.004 

0.001 

- 0.003 

- 0.183 




-0.003 

0.079 

-0.002 

0.012 




-0.022 

o.o46 

-0.004 

0.021 




-0.007 

-0.005 

C.001 

0.222 




0 

0 

0 

0 

0 



0 

0 

0 

0 

0 



- 0.003 

-0.031 

0.001 

0.008 

1.618 


C 

1.064 

0.231 

- 0.061 

-0.031 

- 0.129 



- 0.035 

0.161 

i.n4 

-0.003 

- 0.386 



-0.101 

0.663 

1.443 

-O.O 65 

- 1.562 

35 


0.013 

-1.624 

0.602 

2.218 

0.238 


0 

1.0 

0 

0 





0 

0 

1.0 

0 




0.372 

0 

0.001 

0.007 



D 

-0.011 

-0.002 

- 0.005 

- 0.105 




-0.023 

0.068 

0.002 

0.008 




- 0.081 

0.133 

0.020 

0.041 




0.007 

0.009 

0.121 

- 0.038 




0 

0 

0 

0 

0 



0 

0 

0 

0 

0 



- 0.037 

0.031 

- 0.016 

-0.042 

1.368 

47 

C 

1.081 

0.149 

- 0.057 

0.001 

- 0.086 



- 0.006 

0.098 

1.298 

-0.014 

- 0.236 



-0.053 

0.241 

0.339 

0.076 

- 0.925 



-0.250 

- 1.154 

0.415 

1.751 

0.226 
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'jlABLS t U (Concluded) 


HA-deg l^atrix 


tSatrix Elements 




0 

1.0 

0 

0 




0 

0 

1.0 

0 




0.546 

0.003 

- 0.005 

0.018 


47 

D 

-0.013 

-0.002 

0 

-0.036 




- 0.026 

0.152 

0,004 

0.009 




-0.077 

0.078 

0.018 

0.021 




-0.007 

0.018 

0.108 

0.116 




0 

0 

0 

A 

0 



0 

0 

0 

0 

0 



-0.020 

0.037 

- 0.005 

- 0.025 

1.357 


c 

1.059 

0.119 

- 0.051 

0 

- 0.066 



- 0.016 

0.032 

1.479 

0.025 

- 0.136 



-0.151 

0.124 

0.314 

0.120 

-0.780 

60 


-0.847 

- 0.109 

- 0.126 

1.456 

-0.119 


0 

1.0 

0 

0 





0 

0 

1.0 

0 




0.783 

0.003 

- 0.006 

0.002 



D 

- 0.015 

- 0.003 

-0.004 

-0.064 




- 0.016 

0.192 

0.007 

0.009 




- 0.063 

O.C64 

0.020 

0.064 




-0.047 

-0.012 

0.024 

0.252 




























TABLE Tv 


NORMALIZED POSITION AND RATE 
LIMITS FOR ENGINE CONTROLS 


Iformalized Limits 


Controls 

Rate 

Limits 

Upper 

Position 

Limit 

Lower 

Position 

Limit 

Jet Exhaust Area 

= 1.28 

2.3 

0.8 

Fan Inlet Guide Vanes 

±1.25 

0 

-0.5 

Rear Compressor Variable Vanes 

i 5.00 

0.2 

-2.0 

Main Burner Fuel Flow 

±1.76 

1.2 

0.1 
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TABLE V 

PERFORMANCE INDEX WEIGHTING FACTORS FOR DESIGN OF 
OPTIMAL MULTIVARIABLE ETS3INE CONTROLLER 

Performance Index and Weighting Factors of Eq. ( 3 ) 


Power Lever 
Angle. 
PLA-deg 



Output Weighting Factors 


» 1 

a -2 ^3 j % % 


Control WeiRhting Factors 


r l r 2 


20 

160 

4 o 

ID 

2500 

35 

160 

4 o 

10 

5000 

47 

160 

4 o 

500 

Uoooo 

60 

160 

4 o 

500 

20000 

73 

160 

4 o 

1000 

500 



100 1000 

0 25 


25 16 
25 16 
160 16 


k weighting on perturbed fan inlet guide vanes 
q 2 = weighting on perturbed rear compressor variable vanes 
k weighting on perturbed thrust 
q^ k weighting on perturbed high turbine inlet temperature 
q^ = weighting on rate-of-change of fan stability margin 
qg k weighting on rate-of-change of high compressor stability margin 
r-^ k weighting on rate-of-change of exhaust area 
r 0 k weighting on rate-of-change of fan inlet guide vanes 



r^ = weighting on rate-of-change of rear compressor variable vanes 


k weighting on rate-of-change of fuel flow 


irs ir\ 































TABLE VI 


FEEDBACK GAM H AND L MATRICES FCR 
OPTIMAL MULTIVARIABLE ENGINE CONTROLLER 

Intermediate Values of H and L Scheduled 
Linearly with High Compressor Speed, Ng 


Power 

Lever 

Angle, 

PLA-deg 

*2 

Matrix 

Matrix Elements 

20 

0.6 9 

H 

0.001 0.001 0.073 - 0.057 

- 3.144 - 0.024 0.002 - 1.307 

0.003 - 0.282 -0.006 0.220 

0.064 - 0.057 0.012 - 4.930 

L 

-o.ooi -0.006 0.013 o.o 46 o.on 
- 0.346 -3.457 0.122 1.699 0.055 
- 0.046 -0.494 -0.037 0.907 -0.018 
-1.447 -7.308 0.207 3.635 0.107 

35 

0.83 

H 

-0.027 0.017 -0.481 0.685 

-3.161 0.002 0.019 0.064 

0 -0.627 -0.022 0.283 

-0.037 -0.163 0.124 -31.369 

L 

-0.007 - 0.027 -0.121 -0.452 -0.052 
-o.oo 4 0.047 -0.013 -0.038 - 0.004 
-0.010 0.091 -0.037 -0.094 -0.018 
-0.801 -2.518 1.101 2.756 0.518 

47 

0.88 

H 

0.826 0.505 3.578 - 2.175 

- 2.554 0.086 2.497 0.457 

4.252 - 4.777 8.690 37.490 

0.467 - 0.293 0.223 - 87.033 

L 

-o.ooi 0.209 0.494 1.063 2.302 
0.151 - 0.359 - 0.015 0.752 2.156 
1.254 -2.308 - 2.333 3.453 14.077 
-I.187 - 0.523 - 0.124 0.276 1.860 


(Continued) 
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TABLE VII 

LINEARIZED INLET DYNAMICS 
AT 40,000 FT, MACH 2.2 OPERATING CONDITION 


Matrix 

Matrix Elements 

A 

-100.00 0 

236.40 -37.04 

B 

0 30.00 136.89 

-427.20 0 0 

C 

0 0 

- 0.18 1.00 

0.03 - 0.16 

D 

000 

000 

cvi oj 6 

1 


State Ordering 

Control Ordering 

Output Ordering 

Inlet state 1 (see Fig. ll) 
Inlet state 2 (gee Fig. ll) 

Throat area 

Bypass area 

Fan corrected airflow 

Throat Mach number 

Normal shock position 

Fan inlet total pressure 
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TABLE VIII 

NORMA LIZE D STEADY-STATE INIET-ENGINE PARAMETERS 
AT 40,000 FT, MACH 2.2 OPERATING CONDITION 


Parameter 


Normalized j 

Type 

Parameter Name 

Value 


Thrust 

1.48 


High Turbine Inlet Temperature 

1.07 


Fan Corrected Airflow 

0.71 

Output 

Fan Stability Margin 

0.23 


Compressor Stability Margin 

0.19 


Throat Mach Number 

1.30 


Inlet Normal Shock Position 

0.16 


Fan Inlet Total Pressure 

1.0 


Fan Turbine Inlet Temperature 

1.08 


Main Burner Pressure 

0.88 

State 

Fan Speed 

C.97 


High Compressor Speed 

1.03 


Afterburner Pressure 

0.91 


Jet Exhaust Area 

1.83 


Fan Inlet Guide Vanes 

- 0.50 


Rear Compressor Variable Vanes 

-0.34 

Control 

Main Burner Fuel Flow 

0.85 


Afterburner luel Plow 

1.0 


Throat Area 

1.0 


Bypass Area 

— 

o.ll 


I 
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TABLE 3X 

LINEARIZED engine dynamics 
AT 40,000 FT, MACH 2.2 OPERATING CONDITION 

HA. - 120 deg 


Matrix 

Matrix Elements 


- 55-478 

2.640 

-5.583 

-13.881 

0.651 



27.172 

- 93.200 

48.331 

117.200 

11.393 


A 

0.801 

4.783 

- 5.872 

- 0.351 

- 2.824 



2.146 

1.636 

-0.656 

- 3.967 

-0.560 



3.148 

0.354 

7.076 

2.516 

-8.539 



0.500 

- 0.417 

- 3.085 

28.630 

0 

-21.883 


-0.806 

1.214 

12.911 

4.734 

-0.260 

63.993 

3 

- 0.074 

-O.326 

-0.006 

- 0.275 

-0.045 

- 1.154 


-0.080 

-0.026 

-0.270 

- 0.755 

-0.015 

-0.310 


- 4.362 

0.583 

0.326 

- 1.237 

1.946 

6.195 


0.001 

-0.007 

- 1.840 

- 0.131 

3.405 



1.128 

0.073 

-0.050 

0.064 

-0.058 


c 

o.oil 

- 0.015 

1.168 

0.102 

-0.092 



0.093 

-0.081 

0.658 

0.639 

-0.429 



-0.210 

-0.687 

0.388 

1.451 

0.115 



1.481 

- 0.144 

- 0.017 



-1.270 


-0.009 

0.001 

0.008 

-0.086 

0 

0.061 

D 

-0.007 

0.091 

0.014 



0.108 


- 0.031 

0.078 

0.079 


-0.027 

0.523 


0.004 

0.018 

0.075 

0.063 


0.439 


State Ordering 

Control Ordering 

Output Ordering 

Fan turbine inlet 
temperature 

Main burner pressure 

Fan speed 

High compressor speed 
Afterburner pressure 

Jet exhaust area 

Fan inlet guide vanes 

Rear compressor variable 

vanes 

Main burner fuel flow 
Afterburner fuel flow 

Fan inlet total pressure 

Thrust 

High turbine inlet 
temperature 

Fan corrected airflow 
Fan stability margin 
Compressor stability 
margin 
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TABLE X 

LINEARIZED INLET-ENGINE DYNAMICS 
AT Lo,000 FT, MACH 2.2 OFERATING CONDITICN 


{•Matrix 

Matrix Elements 


-55.478 

-2.640 

-5.583 

-13.881 

0.651 

-0,643 

3.501 


27.172 

-93.200 

48.331 

117.200 

11.393 

1.881 

-10.239 


0.801 

4.783 

-5.872 

-0.351 

-2.824 

-0.034 

0.185 

A 

2.146 

1.636 

-0.656 

-3.967 

-0.560 

-0.009 

0.050 


3.148 

0.354 

7.076 

2.516 

-8.539 

0.182 

-0.991 


1.445 

-2.120 

159.940 

14.024 

-12.610 

-99.565 

-2.370 


0 

0 

0 

0 

0 

236.400 

-37-040 


0.500 

-.417 

-5.085 

28.630 

0 

-3.092 

0 


-0.806 

1.214 

12.911 

4.734 

-0.260 

9.044 

0 

B 

-0.074 

-0.326 

-0.006 

-0.275 

-0.045 

-O.163 

0 

-0.080 

-0.026 

-0.270 

-0.755 

-0.015 

-0.044 

0 


-4.362 

0.583 

0.326 

-1.237 

1.946 

0.875 

0 


-1.024 

12.512 

1.870 

0.501 

-0.338 

2.093 

30.000 


0 

0 

0 

0 

0 

-427.200 

0 


0.001 

-0.007 

-1.840 

-0.131 

3.405 

-0.037 

0.203 


1.128 

0.073 

-0.050 

0.064 

00 
i r\ 
O 

O 

i 

0.002 

-0.010 


0.011 

-0.015 

1.168 

0.102 

-0.092 

0.003 

-0.017 

C 

0.093 

-0.081 

0.658 

0.639 

-0.429 

0.015 

-0.084 


-0.210 

-0.687 

0.388 

1.451 

0.115 

0.013 

-0.070 


0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

-0.184 

1.000 


0 

0 

0 

0 

0 

0.029 

-0.160 


1.481 

-O.l uu 

-0.017 

-o.oo Q 

0.004 

-0.180 

0 


-0.009 

0.001 

0.008 

-0.086 

0 

0.009 

0 


-0.007 

0.091 

0.014 

0.004 

-0.002 

0.015 

0 

D 

-0.031 

0.078 

0.079 

0.001 

-0.027 

0.074 

0 

0.004 

0.018 

0.075 

0.063 

0.006 

0.062 

0 


0 

0 

0 

0 

0 

-2.520 

0 


0 

0 

0 

0 

0 

2.413 

0 


0 

0 

0 

0 

0 

0.l4l 

c 


(Continued) 
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TABLE X (Concluded) 


State Ordering 

Control Ordering 

Output Ordering 

Fan turbine inlet 
temperature 

Main burner pressure 

Fan speed 

High compressor speed 
Afterburner pressure 

Inlet state 1 (see Fig. 11) 
Inlet state 2 (see Fig. 11) 

Jet exhaust area 

Fan inlet guide vanes 

Rear compressor variable 

vanes 

Main burner fuel flew 
Afterburner fuel flow 
Throat area 

Bypass area 

Thrust 

High turbine inlet 
temperature 

Fan corrected airflow 

Fan stability margin 
Compressor stability 
margin 

Throat Mach number 

Normal shock position 

Fan inlet total pressure 
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TABLE XI 

FERFC8MARCE EfDEX WEIGHTING FACTORS FOR 
DESKS? OF INLET-ENGINE COIiTROLLERS 


Weighting 

Factor 

Perturbed Parameter Weighted 

Value 


High Turbine Inlet Temperature 

500 

% 

Fan Corrected Airflow 

6400 

q 3 

Fan Stability Margin 

100 

% 

Compressor Stability Margin 

200 

a -5 

Fan Speed 

500 

% 

High Compressor Speed 

250 

°7 

Throat Mach Number 

100 

% 

Normal Shock Position 

100 

r l 

Rate-of-Change of Exhaust Area 

15 

r 2 

Rate-of-Change of Fan Inlet Guide Vanes 

16 

r 3 

Rate-of-Change of Rear Compressor Variable Vanes 

1 

r U 

Rate-of-Change of Main Burner Fuel Flow 

8 

r 5 

Rate-of-Change of Afterburner Fuel Flow 

10 5 

r 6 

Rate-of-Change of Throat Area 

22 

r 7 

r 8 

Rate-of-Change of Bypass Area 

Fan Corrected Airflow (Inlet Control Design Only) 

2 

10 7 
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Matrix Elements 

-0.185 

-0.151 

- 1.826 

-3.341 

O .832 

-0.012 

-0.007 

-0.137 

-O .189 

-3-946 

-1.240 

0.740 

-0.054 

- 0.012 

-1.844 

- 2.190 

-36.135 

-27.566 

10.438 

- 0.130 

-0.033 

-0.736 

-0.473 

-6.995 

-7.335 

2.566 

-0.094 

-0.044 

0 

0 

0 

0 

0 

0 

0 

0.221 

0.266 

6.092 

1.268 

- 0.907 

0.768 

0.215 

-0.712 

- 0.941 

-19.702 

- 6.188 

3.484 

-1.543 

-0.479 

-1.579 

- 0.093 

-0. 350 

-1.034 

1.021 

0.189 

- 0 . * ’ 

_0.oPP 

_ n *— 

1 


‘ /O 5 

0.659 

- . ‘ • 

-5.245 

- 2.350 

-4.518 

-4.535 

3.541 

2.768 

- 1.397 

-1.939 

- 0.392 

-0.567 

-5.734 

1.257 

1.499 

- 0.301 

0 

0 

0 

0 

0 

0 

0 

0.129 

0.480 

0.126 

0.545 

-0.127 

-13.872 

1.964 

- 0.580 

- 1.565 

-0.699 

- 1.206 

0.537 

21.602 

-6.964 
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TABLE X1IT 

FEEDBACK GAIN MATRICES AND OF THE SEPARATE 
INLET-ENGErE CONTROLLERS FOR WEIGHTINGS SHCWN H? TABLE XI 


Matrix 

Matrix Elements 


- 0.179 

- 0.144 

-1.679 

-3.266 

0.797 

0 



- 0.113 

-0.156 

-3.230 

-1.079 

0.628 

0 



-1.692 

-1.996 

-32.509 

- 26.146 

9.682 

0 


G 1 

- 0.693 

-0.417 

-5.371 

- 6.888 

2.338 

0 



0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0.870 

0.258 


0 

0 

0 

0 

0 

- 1.741 

-0.570 


- 1.575 

-0.082 

-0.340 

- 1.018 

1.018 

0 

0 


- 0.077 

-0.674 

-0.129 

-0.152 

0.080 

0 

0 


- 5.107 

-2.065 

-4.362 

- 4.104 

3.421 

0 

0 

G 2 

-1.908 

-0.304 

-0.513 

-5.680 

1.226 

0 

0 


0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

- 14.984 

2.167 


0 

0 

0 

0 

0 

23.838 

-7.266 


4 o 
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FIG. 3 



IGV. RCVV 























NORMALIZED F401 THRUST RESPONSE FOR OPTIMAL AND CONVENTIONAL CONTROLLERS 



TIME - SEC 








4iu***tm 
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FIG. 5 


I 


I 


I 


NORMALIZED F401 TEMPERATURE AND STABILITY MARGIN RESPONSES 
FOR OPTIMAL AND CONVENTIONAL CONTROLLERS 


RESPONSES PRODUCED BY STEP CHANGE iN PlA FROM 20 DEG (IDLE-9 PERCENT) TO 73 DEG (MILITARY-100 PERCENT) 

OPTIMAL MULTIVARIABLE — — —-CONVENTIONAL CONTROLLER 

CONTROLLER 


(a) HIGH TURBINE INLET TEMPERATURE 



(b) COMPRESSOR STABILITY MARGIN 



(c) FAN STABILITY MARGIN 



N02-89-4 
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N911620-2 FIG. 6 

NORMALIZED F401 SPOOL SPEED AND AIRFLOW RESPONSES 
FOR OPTIMAL AND CONVENTIONAL CONTROLLERS 


I 

I 

I 

I 

I 

I 

I 

I 


RESPONSES PRODUCED BY STEP CHANGE IN PLA FROM 20 DEG (IDLE-9 PERCENT) TO 73 DEG (MILITARY—100 PERCENT) 

— - - OPTIMAL MULTIVARIABLE-—-CONVENTIONAL CONTROLLER 

CONTROLLER 


(a) FAN SPEED 



(b) COMPRESSOR SPEED 



(c) AIRFLOW 



N 02-89-5 







JET EXHAUST AREA 
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FIG. 7 


NORMALIZED F401 EXHAUST AREA RESPONSE 
FOR OPTIMAL AND CONVENTIONAL CONTROLLERS 


EXHAUST AREA RESPONSE PRODUCED BV STEP CHANGE IN PLA FROM 
20 DEG (IDLE—9 PERCENT) TO 73 DEG (MILITARY-100 PERCENT) 


OPTIMAL MULTIVARIABLE — CONVENTIONAL CONTROLLER 

CONTROLLER 
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NORMALIZED F401 FAN INLET GUIDE VANE RESPONSE 
FOR OPTIMAL AND CONVENTIONAL CONTROLLERS 


FAN INLET GUIDE VANE RESPONSE PRODUCED BY STEP CHANGE IN PLA 
FROM 20 DEG (IDLE-9 PERCENT) TO 73 DEG (MILITARY-100 ’ERCENT) 


OPTIMAL MULTIVARIABLE-— — _ CONVENTIONAL CONTROLLER 

CONTROLLER 


INCREASING FAN FLOW AREA 
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FIG. 10 


NORMALIZED F401 REAR COMPRESSOR VARIABLE VANE RESPONSE 
FOR OPTIMAL AND CONVENTIONAL CONTROLLERS 


REAR COMPRESSOR VARIABLE VANE RESPONSE PRODUCED BY STEP CHANGE IN PLA 
FROM 20 DEG (It -9 PERCENT) TO 73 DEG (MILITARY-100 PERCENT) 


OPTIMAL MULTIVARIABLE-CONVENTIONAL CONTROLLER 

CONTROLLER 



HIGH COMPRESSOR SPEED 


N02-89-6 
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FIG. 12(a)-(d) 


NORMALIZED INLET-ENGINE RESPONSE FOR INTEGRATED 
AND SEPARATE MULTIVARIABLE CONTROLLERS 

PERTURBATIONS A80UT NORMALIZED STEADY-STATE VALUES PRODUCED 
BY SIMULATED AFTERBURNER IGNITION 

- RESPONSE FOR INTEGRATED CONTROLLER 

-RESPONSE FOR SEPARATE CONTROLLER 


(a) PERTURBED THROAT MACH NUMBER 



(c) PERTURBED THRUST 



TIME-SEC 


N02-89- 
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FIG. 12 (e)-(i) 


NORMALIZED INLET-ENGINE RESPONSE FOR INTEGRATED 
AND SEPARATE MULTIVARIABLE CONTROLLERS 


PERTURBATIONS ABOUT NORMALIZED STEADY-STATE VALUES 
PRODUCED BY SIMULATED AFTERBURNER IGNITION 

-. — ■ RESPONSE FOR INTEGRATED CONTROLLER 

-RESPONSE FOR SEPARATE CONTROLLER 
























11620-2 


FIG. 12 {JJ—(0 


NORMALIZED INLET-ENGINE RESPONSE FOR INTEGRATED 
AND SEPARATE MULTIVARIABLE CONTROLLERS 


PERTURBATIONS A80UT NORMALIZED STEADY-STATE VALUES 
PRODUCED BY SIMULATED AFTERBURNER IGNITION 


• RESPONSE FOR INTEGRATED CONTROLLER 
RESPONSE FOR SEPARATE CONTROLLER 


\m 


(j) PERTURBED THROAT AREA 


- 0.001 




(k) PERTURBED BYPASS AREA 



(1) PERTURBED JET EXHAUST AREA 



TIME-SEC 


N02-89-8 


CONT'D 









NORMALIZED INLET-ENGINE RESPONSE FOR INTEGRATED 
AND SEPARATE MULTIVARIABLE CONTROLLERS 


PERTURBATIONS ABOUT NORMALIZED STEADY-STATE VALUES 
PRODUCED BY SIMULATED AFTERBURNER IGNITION 

-RESPONSE FOR INTEGRATED CONTROLLER 

-RESPONSE FOR SEPARATE CONTROLLER 


(m) PERTURBED MAIN BURNER FUEL FLOW 





(o) PERTURBED REAR COMPRESSOR VARIABLE VANE POSITION 





















N911620-2 


riu. io 


NORMALIZED INLET-ENGiNE RESPONSE FOR INTEGRATED AND SEPARATE CONTROLLERS 
WITH INCREASED WEIGHTING OF SHOCK POSITION IN INLET PERFORMANCE INDEX 


PERTURBATIONS ABOUT NORMALIZED STEADY-STATE VALUES 
PROOUCEDBY SIMULATED AFTERBURNFR IGNITION 


RESPONSE FOR INTEGRATED CONTROLLER 

RESPONSE FOR SEPARATE CONTROLLER 
WITH INCREASED SHOCK POSITION WEIGHTING 


0.005 

0 

-0.005 

- 0.010 

-0.015 


(a) PERTURBED THROAT MACH NUMBER 


—^-n— 0.0018 


»£*—0.0110 

1 

J_1_ 


(b) PERTURBED NORMAL SHOCK POSITION 



(c) FERTURBED THROAT AREA 



(d) PERTURBED BYPASS AREA 



2.0 

TIME - SEC 


N02-89-9 








